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Abstract. We analyse the constraints obtained from new atomic clock data on the possible 
time variation of the fine structure 'constant' and the electron-proton mass ratio and show 
how they are strengthened when the seasonal variation of Sun's gravitational field at the 
Earth's surface is taken into account. 



General relativity and the standard model 
of particle physics depend on some 27 
seemingly independent numerical parameters. 
These include the fine structure constants with 
determine the strengths of the different forces, 
matrix angles and phases and the relative 
masses of all known fundamental particles. 
Together, these parameters are commonly re- 
ferred to as the 'fundamental constants of 
Nature', although many modern proposals for 
fundamental physics predict that they are nei- 
ther strictly fundamental nor constant. 

Indeed, variation of the traditional 'con- 
stants', at some level, is a common predic- 
tion of most modern proposals for fundamen- 
tal physics beyond the Standard Model. For 
instance, if the true fundamental theory ex- 
ists in more than four space-time dimensions, 
the constants we observe are merely four- 
dimensional 'shadows' of the truly fundamen- 
tal higher dimensional constants. The four di- 
mensional constants will then be seen to vary 
as the extra dimensions change slowly in size. 
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Searches for, and limits on, the variations of 
the traditional constants therefore provide an 
important probe of fundamental physics. 

In many theoretical models of varying con- 
stants one expects the values of the con- 
stants to evolve slowly as the Universe ex- 
pands. Locally, this would manifest itself as 
a slow temporal drift in the values of the 
constants. Laboratory constraints on such a 
drift are generally found by comparing clocks 
based on different atomic frequency standards 
over a period of one or more years. The cur- 
rent best limit on the drift of the electromag- 
net ic fine structure cons tant, a, was found 
by iRosenband et alj d2008l) . The ratio of alu- 
minium and mercury single-ion optical clock 
frequencies, fAi+/fiig+, was repeatedly mea- 
sured over a period of about a year and these 
measurements gave: a/a = (-5.3 ± 7.9) x 
10- 17 yr _1 . 

If the 'constants', such as a, can vary, then 
in addition to a slow temporal drift one would 
also ex pect to see an annual modulation in their 
values (Shaw 120071) . In many theories, the Sun 
perturbs the values of the constants by a factor 
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roughly proportional to the Sun's Newtonian 
gravitational potential, which scales as the in- 
verse of distance, r, between the Earth and the 
Sun. Since r fluctuates annually, reaching a 
minimum at perihelion in early January and a 
maximum at aphelion in July, the values of the 
constants, as measured here on Earth, should 
also oscillate in a similar seasonal manner. 
Moreover, in many models, this seasonal fluc- 
tuation is predic ted to dominate over any linear 
tempo ral drift (IShaw 1 120071: iBarrow & Shaw I 
120081) . 

We suppose that the Sun creates a distance- 
dependent perturbation to the measured value 
of a coupling constant, C, of amplitude 5 In C — 
C(r). If this coupling constant is measured on 
the surface of another body (e.g. the Earth) 
which orbits the first body along an elliptical 
path with semi-major axis a, period T v , and ec- 
centricity e<l, then to leading order in e, the 
annual fluctuation in C, £C an nuai will be given 
by 



annual 



c 



-Cc cos 



2nt 



+ 0(e 2 ), 



(1) 



where cq = eaC'(a), C'(a) = dC(r)/dr\ r=a 
and t = nTp, for any integer n, corresponds 
to the moment of closest approach (perihe- 
lion). For the Earth, a = 149, 597, 887.5 km 
In the case of the Earth moving around the 
Sun, over a period of 6 months from perihe- 
lion to aphelion one would therefore measure 
a change in the constant C equal to 2cc- As 
we noted above, in most theoretical models of 
varying constants, <51nC = C(r) oc AU Q {r) 
where U Q (r) - -GM/r is the Newtonian po- 
tential of the Sun. We introduce sensitivity pa- 
rameters kc defined by 6 In C - kcAU e . Hence: 



GM , n 
c c = ea^r -k c = 1.65 x 10"'%. 



(2) 



Once one specifies a theoretical model of vary- 
ing constants, the kc are determined. Hence 
measuring or bounding the kc constrains the 
underlying theory. One generally expects kc ~ 
0(1) in higher dimensional theories such as 
String Theory, however, as we shall see, obser- 
vat ions typically constrain k c <k 1. 

IRosenband et alj d2008l) fitted a linear drift 
in a to their data finding a/a = (-1.6 ± 
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Fig. 1. Frequency ratio /Ai+//H g + as measured by 
Rosenband et al. (2008). The solid black line shows 
the maximum likelihood fit for a seasonal variation. 



2.3 )x 10~ 17 yr~' . We fitted the expected form of 
any annual fluctuation, Eq.Q] to the measure d 
values of f M+ 1 f Hg jBarrow & Shawl 120081) . 
Fig. Q] shows the best-fit line superimposed 
on the (Rosenb and et al.| [2008) data. It should 
be noted that the magnitudes of the system- 
atic errors for the middle three data points 
were not verified to the same precision as 
they were from the fir st and last data points 
(Rosenban d et al.l[2008l) . We do not expect to 
this have a great effect on the resulting con- 
straint on k a . Using dln(./Ai+/.fag+) = (3.19 + 
0.008)<5a/a, JRosenband et al]l2008l) . a maxi- 
mum likelihood fit to the data gives 



k a = (-5.4 ± 5.1) x 10" 



(3) 



This limit on k a is almost an 2 order of mag- 
nitude improvement on the previo us hmit of 
k a = (2.5 ± 3.1) x 10~ 6 from by iBlatt et all 
(2008). 

The freque ncy shifts measured by 

IRosen band etal.1 d2008h are not sensitive to 
changes in the electron-proton mass ratio: 
fi = m e /m p . However, measurements of 
optical transition frequencies relative to Cs 
are sensitive to both [i and a, and H-maser 
atomic clocks can detect changes in the light 
quark to proton mass ratio: q = m q /m p . We 
define sensitivity parameters kp and k q for n 
and q respectively. We then find that the Yb + 
frequency measurements of iPeik et alJ J2004I) 
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give Jfc ff +0.51^ = (7.1+3.4)xl(T 6 . Combining 
this limit with that on k a and the limits 
k a + .36^ = (-2.1 + 3.2) x 10" 6 dBlatt et all 
2008), k n + 0.17^ = (3.5 ± 6.0) x 10" 7 
dFortier et all 120071) and k a + .13fc q = 
(1 + 17) x 10- 7 lAshbyetal.l (120071) . we find 
(Barrow & Shaw 112008 " 



k M = (3.9 + 3.1) x 10~ 6 , 
k q = (0.1 + 1.4) x 10" 5 . 



(4) 
(5) 



Again these represent improvements on the 
previous limits of k^ - (—1.3 ± 1.7) x 10~ 5 and 
K = ( -1.9 + 2.7) x 10~ 5 found bv lBlatt et all 
(l2008h . 

Seasonal fluctuations occur in varying con- 
stant theories because the constants depend 
on the the vacuum expectation value one or 
more scalar fields, (pi, which interact with 
normal matter. In the solar system the field 
equations of these scalars typically reduce to 
Poisson equations, similar to that satisfied by 
the Newtonian gravitational potential. We la- 
bel the different constants Cj(<f>i) and then: 



30/ SlnCj 



(6) 



We define yS/ ; = dlnCj/d<pj. For small varia- 
tions of the scalar fields and the constants we 
have: 



V 2 In d 



\ r,j 



Slnp 

dlnC, 

1 ) 



4nGp. 



(7) 



The Newtonian potential, O^, obeys: V 2 Ojv = 
4-nGp, and the are defined by In C, = const+ 
kc&N- It follows that for a body orbiting the 
Sun with p = p Q : 



6lnp G 
' ) 



We define A J (Q) 



= 6\np Q /S\nCj; these can 
be calculated from the standard model us- 
ing information about the composition of the 
Sun. Measurements of \\c, from seasonal varia- 
tions therefore determine the model parameters 
Yii PuPij in a model independent fashion. 



The presence of spatial gradients in the 
scalar fields with couple to normal matter also 
results in new or 'fifth' forces. The magnitude 
of the new force towards the Sun on a test body 
with density po, mass mo at a distance r from 
the Sun is given by: 

F 4 = mtfijjAy— = — 2AAo)- (8) 

;' 

Now depends on AL, = 6lnpo/6lnC j 
which depends on the composition of the test 
mass; for a given test mass these can be cal- 
culated using the standard model. This com- 
position dependence violates the universality 
of free-fall and hence the weak equivalence 
principle (WEP). WEP violation searches mea- 
sure the differential acceleration of two differ- 
ently composed test masses towards the Sun 
which is proportional to S;A<j,AA ! ; A/l ! is the 
difference between the V values for the two 
test masses. Hence WEP violation searches in- 
directly probe the kc r The magnitude of any 
composition-dependent fifth force toward the 
Sun is currently constrained to be no stronger 
than about f ew x 10" 13 times than the gravita- 
tional force dSchlamminger et alj2008l) . To ex- 
tract limits on the individual kc i one must ac- 
curately calculate the AA' parameters for differ- 
ent combinations of test masses and have limits 
from different experiments to eliminate the de- 
generacy between the different kc r 

The constraints from WEP tests indi- 
rectly bound k a . Indeed, they often provide 
the tightest constraint s on k„ dShawl l2007t 
IBarrow & Shawl 120081; IDent 112008b . A recent 
and thorough an alysis of the WEP violation 
constraints on k a (|Dent][2008) found: 

\k a \ < (0.23- 1.4) x 10" 8 

with a similar limit on k q . The uncertainty in 
this limit is to due to how one models nuclear 
structure. Despite these uncertainties, it is clear 
that WEP violation constraints from laboratory 
experiments currently provide the strongest, al- 
beit indirect, bounds on k a and the other kc t ■ 

We now consider the sensitivity that would 
be required of an atomic clock to provide 
tighter constraints on k a than those com- 
ing indirectly from limits on WEP violation. 
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Suppose that the ratio of two transition fre- 
quencies, /a//b, can be measured with a sen- 
sitivity o-f, and that S]n(f A /f s ) = S a Sa/a 
(typically S a ~ 0(1), although some tran- 
sitions exhi bit a g reatly increased sensitivity 
(Flambaum 2006)). The sensitivity to changes 
in or is then given by cr a — o-f/S a . By simulat- 
ing data sets, we found that the sensitivity to k a 
is significantly improved if one makes Af m > 12 
measurements per year (at roughly regular in- 
tervals). With N m > 12, by performing a boot- 
strap linear regression with 10 5 re-samplings 
of the simulated data points, we find that the 
sensitivity, cr^ to k a is roughly: 

cr k « 0.69x10' 



^Ny(N m - 1) 

where N y is the number of years for which data 
is taken. The total number of measurements is 
therefore N y N m . Indirect constraints currently 
have a sens itivity no better than cr k = 2.3 xl0~ 9 
(iDent 1 120081) . This would be surpassed by di- 
rect measurements if cr a < 2.5 y/N y (N m - 1) x 
10~ 19 . For example, if measurements are taken 
every 20 days (or so) over a single year (jV m = 
18, N y = 1), we need cr a < 10" 18 . lFlambauml 
(2006) have noted that first excited states in 
the 229 Th nucleus is particularly sensitive to 
changes in a and fi with S a ~ 10 5 , allowing for 
a cr a < 10~ 20 . Such a sensitivity would limit 
k a at the 10" 11 level - two orders of magni- 
tude better than the current WEP violation con- 
straints. 

In summary: we have shown how new lab- 
oratory constraints on possible time variation 
in the fine structure 'constant' and the electron- 
proton mass ratio can yield more sensitive lim- 
its by incorporating the effects of the seasonal 
variation of the Sun's gravitational field at the 
Earth's surface. This seasonal variation is ex- 
pected in all theories which require that the co- 
variant d'Alembertian of any scalar field driv- 
ing variation of a 'constant' is proportional to 
the dominant local source of gravitational po- 
tential. The recent exp erime ntal results from 
dRosenband et al.ll2008l) and dPeik et alJl2004l) 
have reached the sensitivity of the quasar ob- 
servations of varying a and fi made at high 
redshift and we have shown may soon provide 



stronger bounds on varying constants than con- 
ventional ground-based WEP experiments. 
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